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For an exact chiral symmetry which is spontaneously bro- 
ken at zero temperature, we show that at nonzero temper- 
ature, generally pions travel at less than the speed of light. 
This effect first appears at next to leading order in an expan- 
sion about low temperature. When the chiral symmetry is 
approximate we obtain two formulas, like that of Gell-Mann, 
Oakes, and Renner, for the static and dynamic pion masses. 
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Pions are light because they are (almost) Goldstone 
bosons: in QCD, quarks have an (approximate) chiral 
symmetry of SU(2)e x SU(2) r which is spontaneously 
broken to the usual isospin symmetry of SU (2)y by the 
dynamical generation of a quark condensate Era]. No- 
tably, the pion mass squared is proportional to the up and 
down quark masses through the formula of Gell-Mann, 
Oakes, and Renner 

In this paper we consider how pions propagate in a 
thermal bath ||; similar results should also hold for pi- 
ons propagating in a Fermi sea of nucleons. In the limit 
of exact chiral symmetry, pions are true Goldstone modes 
and so massless. At zero temperature relativistic invari- 
ance then requires pions to travel at the speed of light. 
Our basic point is elementary: since the presence of a 
medium provides a privileged rest frame, relativistic in- 
variance no longer applies, and so typically pions travel 
at less than the speed of light. We also derive how the 
formula of Gell-Mann, Oakes, and Renner generalizes to 
nonzero temperature Because the pion's velocity is 
less than c, in a medium the pion dispersion relation, as 
a function of momentum, is "flattened" from that at zero 
temperature. Such a flattening has been found in a wide 
variety of models [p|-|l0t , due apparently to the detailed 
dynamics. Our results show that at least some of the flat- 
tening arises on very general grounds, as a consequence 
of chiral symmetry breaking in a medium. 

While we speak of pions throughout, our conclusions 
apply to Goldstone bosons in any system which is rel- 
ativistically invariant at zero temperature. Indeed, our 
results for the changes in the pion dispersion relation 
have exact analogies with spin waves in antiferromagnets 

We think that 

our manner of derivation — in 
terms of the pion decay constants — is novel and illumi- 
nating. The detailed calculations which we perform to 
demonstrate this effect in the linear sigma model extend 
and complement previous results by Itoyama and Mueller 

We begin with a heuristic derivation, in the limit of 



exact chiral symmetry. At zero temperature, the matrix 
element of the axial vector current, A£, sandwiched be- 
tween the vacuum and a pion of momentum P M = (p°,p) 
is: 

(Q\AZ\ir b (P))=if w 6 ab P» , (1) 

with a and b isospin indices. The pion decay constant 
f w ~ 93MeV; whenever we write / T , we mean its value 
at zero temperature. 

At nonzero temperature, because of the presence of the 
medium, we expect that there are two distinct pion decay 
constants, one for the timelike component of the current, 

J 7T ' 

(0\Al\K b (P)) T = ifl5 ab p° , 

and one for the spatial, /£, 

<0|A> h (P)) T = z/^V ■ (2) 

Both matrix elements are computed at a temperature T 
in the imaginary time formalism. Implicitly the timelike 
component of the momentum, p , is analytically con- 
tinued from euclidean values (pions, as bosonic fields, 
have p° — 2-nnT for integral n) to Minkowski values, 
p° = —iui + + . In (||), f* and are defined about zero 
momentum, uj and p — > 0. 

The possibility of two distinct pion decay constants is 
familiar from nonrelativistic systems, such as discussed 
by Leutwyler fll2| ; in this context it was recognized previ- 
ously by Kirchbach and Riska and by Thorsson and 
Wirzba §. 

By assumption the chiral symmetry is exact, and only 
broken spontaneously by the vacuum. Consequently, 
while the axial current acts nontrivially on the vacuum, 
it nevertheless is conserved on the pion mass shell. At 
zero temperature this is trivial: the divergence of the ma- 
trix element in ffl) is (0|9 M ^4 M |7r) ~ firP 2 , which vanishes 
when P 2 = —uj +p 2 = 0, as expected for a massless, 
relativistically invariant field. 

At nonzero temperature, however, the condition that 
the axial current is conserved on the pion mass shell 
leads to interesting restrictions on the two pion decay 
constants. The divergence of the axial current in (|^) 
vanishes when 

flPa + fvP 2 = Ok mass shell (3) 

At nonzero temperature each pion decay constant, f* 
and /*, has a real and an imaginary part. The pion mass 
shell then lies in the complex plane, at p° = —iuj — 7. 
Equating the real parts of (||) gives 
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The requirement that pions travel at less than (or equal 
to) the speed of light, v < 1, implies Re/^ < Re/*. To 
obtain this, we assume that the imaginary parts can be 
neglected relative to the real parts, 



Im/*' s «Re/* 
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Physically, v < 1 is most familiar: pions move through 
a medium as if it has an index of refraction greater than 
or equal to one. 

The imaginary part of the mass shell is given by 
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x , /*(T) and f%(T) -> |5|, and the region in which 

About T x , over 



The requirement that pions are damped, and not anti 
damped, fixes 7 to be semi positive definite; using (Q), 
this then constrains the real and imaginary parts of /* 
and 

Our analysis only applies to "cool" pions, where the 
components of the pion momenta, uj and p, are small 
relative to the real parts of j\ and /* . If the chiral phase 
transition is of second order at T — T x , then as T — > 

cool pions dominate shrinks to zero 
large distances the behavior of pions (and the a meson 
is controlled by an 0(4) critical point, as appropriate for 
two massless flavors. 

Assuming that the imaginary parts of /* and /| are 
nonzero at zero momentum, from (^) the damping rate 
vanishes linearly about zero momentum, 7 ~ p as p — > 
0. This is consistent with Goldstone's theorem M: at 
zero momentum, the complete inverse pion propagator 
must vanish, including both the real and the imaginary 
parts. If 7 ~ p, then the imaginary part of the pion 
self energy, ImIi(P) in (10), and so the complete inverse 
pion propagator A _1 (P), vanishes ~ p 2 as p — > 0. This 
implies that even when pions are damped, about zero 
momentum they still dominate the correlation functions 
of axial vector currents. 

In a nonlinear sigma model, the first contribution to 
the damping rate appears at two loop order [R-18|. Us- 
ing a virial expansion, such as (2.4) of ref. ||17|| 7~we es- 
timate that about zero momentum in the chiral limit, 
7 ~ p(T 4 //*). In the linear sigma model considered be- 
low, the damping rate vanishes exponentially, (|2^), but 
this is special to the kinematics at one loop order in this 
model §J]. 

To make our conclusions rigorous, and to extend them 
to an approximate chiral symmetry, we follow Shore and 
Veneziano jj) by using a chiral Ward identity of QCD. 
Take two flavors of quarks, each with a (current) quark 
mass = m. A chiral Ward identity between the form fac- 
tors and the propagators of the quark composite operator 
^=^7 S gis| 

d^\A^\) T + (qq) T (0|T*^^|0) T 1 = 2mS ab , (7) 



where (qq)r is the quark condensate and {O|7"*0g0| |0)^ 1 
the inverse propagator for the (j>§ field. As usual, this 
chiral Ward identity has the same structure as at zero 
temperature, except that now thermal expectation values 
enter. Assume that ^5 is directly proportional to the pion 
field, 

= m . (8) 

The normalization constant "6" is a function of both tem- 
perature and momentum. The temperature dependence 
follows from our analysis, while we neglect any momen- 
tum dependence. As discussed by Shore and Veneziano 
H , dropping this momentum dependence is equivalent to 
the usual assumptions which give the partial conservation 
of the axial vector current. 

Using (g) and (§) in @, the chiral Ward identity be- 
comes 

- b {flpl + ftf) + b 2 {qq) T A- l (P) = 2m . (9) 

For the inverse pion propagator A~ 1 (P) we take 

A-\P) =p 2 Q + v 2 p 2 +ml-ilmU(P) . (10) 

Similar forms of the pion propagator have appeared pre- 
viously [p[-^o|; for us this form is motivated by the need 
to satisfy the chiral Ward identity. 

Ref. P] requires that the pion field is canonically nor- 
malized, so the coefficient of Pq in the pion propagator 
must be unity. We allow for a pion velocity which is 
less than one by introducing the velocity 'V . Since the 
quark mass m / 0, we introduce a pion mass, m n . Lastly, 
we introduce an imaginary part of the pion self energy, 
ImII(P), which is a function of momentum. This form 
of the propogator should be valid in an expansion about 
zero momentum. 

The chiral Ward identity shows that the assumption 
used to derive (^|) is correct: in the chiral limit, m — 
0, the divergence of the axial current vanishes on the 
pion mass shell, as defined by the condition A~ 1 (P) = 0. 
Since the chiral Ward identity holds for arbitrary (small) 
momentum, however, we can derive several identities by 
matching the coefficients of Pq, p 2 , and 1, for both the 
real and imaginary parts. 

Equating the terms ~ p\ fixes the constant of propor- 
tionality between the quark operator and the pion field 
to be 
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Since both terms on the right hand side of (|llj) change 
with temperature, so does the factor "&" . Matching the 
terms ~ p 2 fixes the velocity as in (f|) . Lastly, matching 
the imaginary parts in gives ImII(P) = 2^7, with 7 
as in @. 

Away from the chiral limit, we match the real parts at 
zero momentum, po = p = 0, to obtain the generaliza- 
tion of the relation of Gell-Mann, Oakes, and Renner to 
nonzero temperature: 
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This is the same expression as Dashen || found at zero 
temperature, except that instead of f n , at nonzero tem- 
perature the real part of /£ enters. A relation like jl^ ) 
was obtained by Thorsson and Wirzba 0); they did not 
recognize, however, that in general has an imaginary 
part, and so wrote just instead of Re/^. 

The pion mass in (|l2|) is the dynamic pion mass, de- 
fined as the position of the singularity in the pion prop- 
agator in the complex po plane at p = 0. Alternately, 
we can introduce the static pion mass, as the position of 
the singularity in the pion propagator for po = in the 
complex p plane. From the form of the pion propagator, 



m^/v, and so by (f|) and (|l2] ) this is just 



{m s J atic f = 
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Obviously, v < 1 implies that 



> m„ 



(13) 
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We now consider where these effects first appear in 
an expansion about zero temperature. Using either a 
nonlinear |Eo| or a linear |^l| sigma model, to leading 
order in Tyf%, 



/<(T) = £(T) = 1 - 



T 2 
12/* 



(15) 



Hence to leading order in low temperature, pions move 
at the speed of light and are undamped. This was es- 
tablished by Dey, Eletsky, and Ioffe p2| , who showed 
that to - T 2 // 2 , the thermal average of the two point 
function of either vector or axial vector currents is di- 
rectly proportional to a linear combination of those at 
zero temperature. Since these two point functions are 
lorentz covariant at zero temperature, they remain so to 
~T 2 // 2 . 

Thus the first place where the effects which we are 
discussing can enter is at next to leading order, ~ T 4 
|p3| . The pion damping rate |l(J and self energy [ p7|Jl8| 
have been computed to ~ T yf% in a nonlinear sigma 
model . In particular, Schenk (l7]]l^] computed the pion 
self energy not in the chiral limit, but using physically 
reasonable approximations. His results imply that for 
T = 150MeV, v ~ .87 @. 

Instead of computing to two loop order in a nonlinear 
sigma model, to illustrate the effect we calculate, in weak 
coupling, to one loop order in a linear sigma model. In 
cuclidean space time the lagrangian is 



— ha 



(16) 



where <j) = (f, 7?) is an 0(A) isovector field. We introduce 
a background magnetic field h which is proportional to 



the current quark mass m. For h — 0, the vacuum ex- 
pectation value of the a is cto = \J p, 2 / A, where we then 
shift a — > (To + c; for two flavors, /„• — a^. 

We compute terms of order ~ T 4 /(/ 2 m 2 ). There are 
many terms of order ~ T 4 // 4 , at both one and two loop 
order. In weak coupling, however, the a meson is light 



relative to / ff 



= 2A/ 2 . Thus in weak coupling, which 



we assume, the terms of order ~ T 4 // 4 are smaller by 
~ A than those computed. Conversely, in the limit of 
strong coupling, m a — > oo, the only terms are those ~ 
TV ft- , 

The diagrams for the pion self energy have been com- 
puted to - T 4 /(/ 2 m 2 ) by Itoyama and Mueller 
/* = /* has been computed to ~ T 2 // 2 by Bochkarev 
and Kapusta l^lj , so all we have to do is extend the cal- 
culation of /* and /| to this order. Consequently, we 
merely sketch the simplest way of performing the calcu- 
lations. For the real parts, the terms of interest arise 
from diagrams involving a virtual a and a tt in a loop, 
such as 



1(P) = tr K 



1 



K 2 ((P - K) 2 +m 2 ) 



(17) 



where tr K =TJ^_ 00 J d 3 fc/(27r) 3 . To - T 2 , it suffices 
to approximate this integral by its value at zero momen- 
tum, neglecting the K dependence in the a propagator, 
so ( |l7] ) becomes 



J(P) 



t TK lo 



ri 2 12 



(18) 



In the integral we have ignored apparent ultraviolet di- 
vergences to concentrate on the term ~ T 2 . Of course 
renormalization is taken care of as usual at zero temper- 
ature. 

To compute terms of ~ T 4 , it is necessary to expand 
the integral in ( |l7j ) to ~ P 2 , including both terms ~ P 2 
and terms ~ P^P V . Besides the integral in (|l8|), we also 
need 



tr K - 



(5^ - ^n v ) 



n 2 T 4 
90 



(19) 



where n M = (1, 0). 

The imaginary part of expressions cannot be extracted 
so easily. We evaluate the imaginary part only near the 
pion mass shell, which is for u> ~ p. In this region, the 
only contribution to the imaginary part of (181) is from 



Im I(P) = 



d 3 k Tt(rii — n 2 ) 
(2tt) 3 4Ei£ 2 



5 (u + Ei - E») . (20) 



In this expression E\ = k is the energy of the pion, E% = 
yip — k) 2 + m 2 is the energy of the a, and n\ = n(Ei), 
n 2 = n {E2) are the corresponding Bose-Einstein distri- 
bution functions. This result can be obtained in various 
ways, such as following p5| j. In all there are four possi- 
ble S- functions in energy which contribute to Im T{P). 



3 



For uj ~ p <C T only that in ( |20[ ) contributes, and corre- 
sponds to Landau damping. In this region, the <5-function 
requires 



(21) 



2(u+pcos9) 



We assume that k 3> m a , and then expand the ener- 
gies accordingly; this is justified, since from (^l|), when 
m a ^> u>,p, then /c 3> m a . The result for the imaginary 
part is 



Im I(P)U~ p<m „ ~ exp {-^ T j 



(22) 



Because the fields being scattered have large momentum, 
the Bose-Einstein distribution functions are essentially 
Boltzman, which generates the exponential suppression 
seen in (|22|) . 

These integrals are sufficient to reproduce the results 
of ref. jl3) for the pion self energy. To evaluate the corre- 
sponding terms for the pion structure constants, we need 
the axial current in the linear sigma model, 



A a = (<to + o)d^ a - 7r a d„a 



(23) 



The diagrams which contribute at one loop order to /* 
and /* are given in fig. (5) of [ pl| . Besides the pion 
self energy, there is a contribution from a ct-tt loop at 
the vertex for A% . These contributions can be evaluated 
expanding integrals like ( |l7|) and using ([To|) . For the 
imaginary parts, we need the integrals 



tr K 



2 \ / ™2 



(24) 



and 



tr K 



k' 



K 2 ((P-K)2 + ml)^ p<<m ° 



16pn \4p J \ ApT J 



(25) 



The results of the computations are as follows. At one 
loop order the quantity 



ti 



T 2 



(26) 



typically enters. To the order we work, we also need 

vr 2 T 4 



t? 



1 mz 



45 flml ' 



327r/^— V IpTj 



exp 



(27) 



Then at weak coupling in the linear sigma model, to 



ft ~ (1 - h + zt 2 + %t 3 ) u 
f°~(l-h-5t 2 -it 3 )U. 
By (Ji|) the pion velocity is 

v 2 - 1 - 8t 2 , 
while from (|^) the pion mass shell is 
ip° ~ vp — ipts . 



(28) 



(29) 



(30) 



Including the one loop self energy computed to this order, 
the pion propagator is 

z^-\p)\^ p ^ ma ~ (i + ti + m 2 ) P 2 + (i + h - 2t 2 ) P 2 

+ ml(l + M 1 /2)-2ip 2 t 3 . (31) 

To ensure that A~ 1 (P) has canonical normalization we 
introduce a factor for wave function renormalization of 
the pion, 



Z n ~ 1 + ti + 6< 2 



(32) 



It is elementary to check that the zero of ([n]) agrees 
with (30). We have included the results to leading order 
in the external field h, when the pion mass is nonzero. 
Assuming that the quark condensate is proportional to 
the vacuum expectation value of the a field, 



(g (Z ) T - ( To(T 1 )- ( To(0)(l-3ti/2) , 



(33) 



we also verify our generalization of the formula of Gell- 
Mann, Oakes, and Renner in ( |l2| ) for the dynamic pion 
mass: 



m2(T)~m£(l + ti/2-6t 2 ) 



(34) 



We conclude with some general comments. First, while 
the effects computed at low temperature ( p(i| ) - (|34|) are 
small, that does not mean that they remain so for tem- 
peratures of physical interest, as seen in the results of p fl - 
JlfH . Secondly, the form of the inverse propagator in ( Jicj ) 
applies not just to Goldstone bosons, but to any scalar 
field at nonzero temperature. For example, in numerical 
simulations on the lattice in euclidean spacetime, typi- 
cally what is measured is only the static mass, not the 
dynamic. 

Finally, we note that the coefficient of v 2 — 1 in (|2£ 
is proportional to the free energy density for pions, 
= 7r 2 T 4 /30. (It would be interesting to know what the 
analogous coefficient is for the nonlinear sigma model in 
the chiral limit.) This and other examples |26| hint of a 
general relation, valid for all temperatures, where the de- 
viation of the velocity squared from unity is proportional 
to the free energy density M. 
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